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Perturbative Coulomb gauge Yang-Mills theory within the first order formalism is considered. 
Using a differential equation technique and dimensional regularization, analytic results for both the 
ultraviolet divergent and finite parts of the two-point functions at one-loop order are derived. It is 
shown how the non-ultraviolet divergent parts of the results are finite at spacelike momenta with 
kinematical singularities on the light-cone and subsequent branch cuts extending into the timelike 
region. 
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^ ' 1. INTRODUCTION 

\ Coulomb gauge QCD is rather special. Amongst all the various gauges, it can be shown in Coulomb gauge that 
the number of dynamical variables reduces to the number of physical degrees of freedom [ij. This allows for a 
I tantalizing glimpse of possible nonperturbative descriptions of confinement and the hadron spectrum. The so-called 
. "Gribov-Zwanziger" scenario of confinement, P, 0]) becomes especially relevant in Coulomb gauge. In this picture, 
(~| ' the temporal component of the gluon propagator provides for a long-range confining force whilst the transverse spatial 

components are suppressed in the infrared. Various calculations support this picture, among them [3, 0J&S[3• 
^^' Perhaps a touch ironically, perturbation theory is one of the starting ingredients for nonperturbative calculations 
D \ in the sense that in the high energy region (and for asymptotically free theories such as QCD this is the perturbative 
1-^^' domain), the regularization and renormalization of the theory play an important role in constraining the necessary 
\ approximations. In^Coulomb gauge, only the leading divergence structure is known due to severe technical difficulties 
T— I ■ (see for example (sj, (9{). For Coulomb gauge within the (standard) second order formalism there exist so-called 
^ ■ energy divergences, which have been shown to cancel up to two-loops @, but a general proof of this cancellation 
Op , is sadly lacking. One way to circumvent the energy divergences is to work within the first order formalism, where 
^ ' formal arguments show that such divergences cancel exactly T|. This circumvention comes at a price: the Dyson- 
, Schwinger equations become cumbersome (lo| and full multiplicative renormalizability is not maintained Whilst 
• ' the leading divergences do provide crucial information about the renormalization of the theory, the remaining finite 
parts are critical to further progress in the field. Having expounded the necessity for perturbative results within the 
nonperturbative context, obviously perturbative results for physical high energy processes are desirable - not the least 
in order to compare with results from covariant gauges. 

The technical barrier to progress in Coulomb gauge perturbation theory stems from noncovariant loop integrals of 
the type: 



o 
o 



X 



/ — (1.1) 

where (in Euclidean space) uP' = uj1+uj'^. Standard techniques such as Schwinger parametrization fail due to the 
complexity of the resulting parametric integrals. One might imagine that using contour integration to firstly evaluate 
the temporal component of the integral might make the situation simpler. However, in such a method, translational 
invariance is lost and since the subsequent spatial integral is ultraviolet [UV] divergent, the result will in general be 
incorrect. A UV-cutoff procedure will also fail. There are however techriiques that can overcome these difficulties and 
one of these is the differential equation technique. In its original form ITsj. complicated massive integrals arising 
in covariant gauge calculations can be considered and when supplemented with integration by parts identities [14„ 1^ , 
the technique becomes a powerful tool. The ethos of the technique is that whilst the multi-dimensional parametric 
form of the integral may be practically impossible to work with, the original integral is itself only a function of a few 
variables and where differential equations can be derived, finding the solution involves integration over only these few 
variables and sorting out the boundary conditions. 

In this paper, we consider the one-loop perturbative two-point functions of Coulomb gauge Yang-Mills theory within 
the first order formalism. Within this noncovariant setting, a variant of the differential equation technique and the 
integration by parts identities are derived in order to evaluate integrals such as the one above. This allows for a full 
analysis of the various propagator and two-point proper functions of the theory. 
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The paper is organized as follows. We start by briefly reviewing the first order formalism used and express the 
two-point functions in terms of their loop integrals. In Section 3, the noncovariant integrals inherent to Coulomb 
gauge are evaluated. This section comprises the bulk of the development necessary to the study and is unashamedly 
technical in nature. The results for the two-point functions are collected in Section 4. We finish with a summary 
and outlook. Those loop integrals that can be evaluated using standard techniques are described in Appendix [XI 
Appendix [Bl contains a nontrivial check on the noncovariant integrals. 



2. THE FIRST ORDER FORMALISM AND PERTURBATION THEORY 



Since Coulomb gauge Yang-Mills theory within the first order formalism is rather different to Yang-Mills theory in 
linear covariant gauges, let us begin by reviewing those aspects of the formalism that will be relevant. For a complete 
description, the reader is referred to Ref. The generating functional is written (in Minkowski space) 



Z[J] = y exp {iSb + iSfp + iS^ + iS,) 



(2.1) 



where $ denotes the collection of fields and the terms in the action are given by 

1 



Sb = 




Sfp = 


j d'^x 




j d^x 



2 



(2.2) 



with the source term defined in condensed notation as (Greek indices such as a refer to all attributes of the field, 
including its type, and summation over all discrete indices and integration over all continuous arguments is implicitly 
understood): 



(2.3) 



In the above, A and a are the spatial and temporal components of the gauge field, tt and arise in the construction 
of the first order formalism (they represent the transverse and longitudinal components of the conjugate momentum 
to the gauge field), c and c are the Grassmann- valued Faddeev-Popov ghost fields introduced by fixing the gauge, A 
and r are Lagrange multiplier fields. The chromomagnetic field, B, is given by 



Bf = e 



ijk 



1 



;9f 



abc j^b j^c 
j k 



(2.4) 



(roman subscripts indicate spatial indices) and the spatial component of the covariant derivative in the adjoint 
representation is 



(2.5) 



The general forms of the various Green's functions that we will be considering are constrained in several ways. 
Since the derivation is necessarily somewhat longwinded, for brevity we omit the details here and again refer the 
reader to Ref. 10] for a full account. There are three constraints. Firstly, the Lagrange multiplier (A, r) and 4> 
field equations of motion can be solved exactly. The former primarily supply the transversality properties of the 
vector-vector propagators (i.e., the connected two-point Green's functions), the latter relating the proper two-point 
functions involving functional derivatives of the i^-field to contractions of those involving the corresponding 7r-field. 
Secondly, the equation stemming from the invariance of the generating functional under the BRS transform (the 
Ward-Takahashi identity in raw functional form) tells us that the A~A propagator must vanish. Thirdly, the discrete 
symmetries of time-reversal and parity constrain the allowed forms with two main consequences: most of the scalar- 
vector propagators must vanish (this is applied in conjunction with the transversality conditions arising from enforcing 
the Lagrange multiplier equations of motion) and the dressing functions of the propagator or two-point proper Green's 
functions must be functions of the variables /cq ^'^d fc^. 

In TablelH the general decomposition of the propagators (collectively denoted by W) is presented. The vector-vector 
propagators are explicitly transverse, with the transverse projector in momentum space given by tij{k) = Sij — kikj / k'^ . 
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It is understood that the denominator factors involving both temporal and spatial components implicitly carry the 
appropriate Feynman prescription, i.e., 



such that the integral over the temporal component can be analytically continued to Euclidean space (fco — > iki). 
Supplemented by the additional expression for the ghost propagator, 

WfAk) - --5"^'^, (2.7) 

the list is complete. Each of the dressing functions Dap is a function of fcg and except the ghost which is a function 
of fc^ only. The tree-level propagators are given by 

Daa = Da. = D,, = D,, = = D„x = D, = 1, D^^ = D^x = 0. (2.8) 

The general decomposition of the proper two-point functions (collectively denoted by F) is given in Tabic |TT1 The 
vector-vector functions contain longitudinal components and the longitudinal projector is written lij{k) — kikj/k"^. 
The ghost proper two-point function is: 

Vf^{k) = S'^^iTS^. (2.9) 

Again, the dressing functions are functions of fcp and fc^ except that for the ghost which is a function of k^ only. At 
tree-level 

TaA = ^AiT ~ Ttttt = rTTCT = Tc = 1, TaA — Tatt = ^TT-n = r^cr = Fo-o- = 0. (2-10) 

The two sets of functions (propagator and two-point proper Green's functions) are related via the Legendre transform 
and we have 

Daa 







Daa 


1 








F2 - 




(fe^ 


-fc^)rAA 


Drjitt, 




r 






p2 




(fc^ 


-P)rA„ 


Drj<t, 




V 




-fc^r^Ar,^)' 


F2 - 

TTtT 




DcTc -- 


= 1, 









= Daa^Acr — D^^ (F^tt + ^ A-k) + Da\, 
= —Da4,TAa — Dcj,^ {Tatt + ^ Air) + D^\, 
= Faa - § [DaxTAa + D^x [V A. + T A.)] ■ (2.11) 

In addition to the two-point functions we have various three-point functions. The tree-level vertices (three-point 
proper Green's functions) used in this work are given by (all momenta are defined as incoming): 

p(0)a6c _ _„fabcc 
TTcrAij — yJ "ij' 

^3AijkiPaiPb,Pc) = -Wf^" [^y (Pa - Pb)k + 5.jkij>b - Pc)t + 5ki{Pc " Pa)j] , 

riSf = [r^f''- ~ fadeJcbe^^ ^ g^^g^^ ^fabejcde „ jade fbce^ ^ g^^g^^ ^jacejdbejabejcde^ | ^ 

4°!r^(p^'Pc,PA) = -^gr'''p-c^, 

T^:}tip,,P.,PA) = ^P,,T^Ta% - ^wr'^P,.. (2.12) 

The (Minkowski space) loop integration measure for the integrals entering the Dyson-Schwinger equations is 
(— dujM ) where 



" (2^)-^+! ^^-^^^ 
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f .(h.\ '•Da A 


f..(u\(-k°)DA^ 













I'liyi^j (i,-2t.2\ 








f (f.\ k"DA^ 
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kj 
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TABLE I: General form of propagators in momentum space. The global color factor 5°'^ has been extracted. All unknown 
functions are dimensionless, scalar functions of and P. 
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TABLE II: General form of the proper two-point functions in momentum space. The global color factor 5°''' has been extracted. 
All unknown functions Fag are dimensionless, scalar functions of feo and P. 



and d = 3 — 2e is the spatial dimension. In order to preserve the dimension of the dressing functions (or, more 
formally, the action) we must also assign a dimension to the coupling through the replacement 

/ ^ 9^^^' (2.14) 

where /U is the square of some non-vanishing mass scale (and which will be later identified with the renormalization 
scale). 

Although the formalism as presented thus far is written in Minkowski space, in order to evaluate the resulting loop 
integrals it is necessary to analytically continue to Euclidean space {ko ik^) and we use the notation to denote 
the temporal component of the Euclidean 4- momentum such that fc^ = fcf + P. The Euclidean integration measure 
is written 

- (^- (2-15) 

Since, in this noncovariant setting, the evaluated dressing functions will be functions of the two variables fc| and fc^, 
the Minkowski space dressing functions are recovered by analytically continuing back with —kQ. Assuming that 

the loop integrals can be expressed in terms of known analytic functions, the Minkowski space dressing functions can 
be given by simply extending the argument fc| to negative values and observing contributions generated by continuing 
through singularities. Physically, such singularities can only occur for lightlike momenta (assuming that there are no 
timelike resonance states) just as in linear covariant gauges. We will use the same notation for the Minkowski and 
Euclidean space functions, since it is clear from the argument or fcl in which space they reside. 

Now, the stated purpose of this paper is to derive the one-loop pcrturbative form for the various two-point dressing 
functions. In linear covariant gauges, the Feynman rules can be easily applied to give directly the expressions for 
propagator functions. The reason for the relative simplicity is that the proper two-point functions are related directly 
to the corresponding propagators via inversion. This is not so within the first order formalism - since the proper 
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two-point Green's functions and the propagators are related effectively by a matrix inversion, one must expand the 
full set of Dyson-Schwinger equations for the proper functions in powers of the coupling, g, and then use the relations 
p. lip to construct the propagators. [Actually, a similar situation exists for the quark propagator in covariant gauges 
too, but the matrix inversion only involves different contractions of the same equation.] In the remainder of this 
section, we will thus consider the one-loop perturbative expansion of the various Dyson-Schwinger equations and 
write the dressing functions of the proper two-point functions in terms of the loop integrals (these integrals will be 
evaluated in the course of the next section). There are seven such equations (derived in fl^, with the exception of 
T^A which is an obvious extension of T^a ) and they read: 



r^^i/c(fc) 
^a-'Amik) 



-^efd 



S^'h -J{- Alum )rTA%ik, k)W'„l{iu)Tlii,{uj, k - cu, -k)wi\^^{iu - fc), 

^5^%k - j{- duJM )rl"jl^ -(fc, ^ - k)W'ff,{u)Tl^^^^,{oj, k - CO, -fc)W^„%(^ - fc), 
- / (- du^M ) r^^hti^ ~ k, k, ~u.)WX%,{u.)Tli'^,,^{^, k - ^, -fc)W^i;,,(a; - k) 



(- dcvM ) ryjX/(^ - fc, k, -u;)WX%,,iu;)rli%,Jcj, k - u;, -k)Wil{u; - fc). 



efd 



ntik) = - / {-du:M)Tff^^,{^^k,k,~u:)W%^,{u:)Tlfl,^^^ 



(- du:M ) ryjr/(^ - fc, k, -u:)Wf,,^{u:)T^^Z^{u:, k - -k)Wi;{u^ - fc). 



efd 



(2.16) 
(2.17) 
(2.18) 
(2.19) 

(2.20) 

(2.21) 



^ AAim 



(fc) 



15'^^ \kH,^ - hkj^ +J{~ diUM )y'i^M%^ - fc, ~io, k)Wf{uo)Tg\^{u, k-uj, -k)Wi\cu - k) 

(- )r2A?('- - k, -u, k)Wfp{u)T'f:^^{u, k - ^k)Wil{uo - k) 
- I {- du^M )rTMA^ k, k)W^j{u;)Tf:^,^ico, k - c., ~k)wit,,{u^ ~ k) 
^ / (- du^M - k, -u;, k)W^%,{u:)Tf:^,^^{u;, k - -k)Wi\^,{u - k) 

(- dcoM)i- dvM )ri°iffct(-«' ^LO,V + CO-k, k)W'^Jk - V - i0)W^A',koi^WZpiv) X 



1 

~2 
1 
"6 

^fghe 
X-ySAnopm 



{k — uj — v,uj, V, —k) 



+ \ I {-du;M)rf}Zfkik,-k,u;,-co)Wl%{-u) 



+- 



(- dcoM){~ dvM )riXt(-«' -uj,v + uj-k, k)WZMWZkoi^n',Lopiv, ~v-u)x 



df 



^fgh 



W^a';,,(« + ^K^A^rmiv + u; , k - V ~ u; , ~k)W^i.^{cj + v-k). 



(2.22) 



These equations are all expanded to one-loop, the contributing terms tabulated in Table IIIII Introducing some 
notation for the tree-level and one-loop terms, we write 



- ap ■ 



2t-(1) 

/3 



(2.23) 



where the dimensionful parameter /i^ is included in F^^^. After inserting the appropriate tree-level factors, resolving 
the color algebra and analytically continuing to Euclidean space, the one-loop expressions read: 



d^jT'il{klk')+k,ik)T'Ai{kie) = -TVe/i^ 

r(i)(fc2) = -N,^^ 



duj LU4 



k^uj'^ik — cj)2 

duj 
uj-^{k^ujf 



LU^UJk 



1 



k-LU^ 



(2.24) 
(2.25) 
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equation 


integral term 


contribution(s) {a, P, ■ . .) 


1\a. (I2.16f 


1st 


(a = 7r,/3 = cr) 


r.. (|2.17P 


1st 


ia = A) 


r^„. (|2.18|) 


1st 




IV.. (|2.19f 


1st 


{a = A,f3 = a) 


IVa, P.20P 


1st 


ia = A,f3 = A) 




2nd 


(q = cr,/3 = tt) 


r,,,,, p.2ip 


1st 


(a = = 7r),(a ^ tt,I3 = A) 




2nd 


— 




1st 


(explicit) 




2nd 


(a = cr, /3 = <;/>) 




3rd 


(a = TT,/? = cr) 


Taa, (|2.22| 


4th 


(q = ^,/3 = A) 




5th 






6th 






7th 





TABLE III: Contributing integral terms to the one-loop perturbative expressions for the Dyson-Schwinger equations (|2.16f - 



a;2(fc-a;)2 



6ik 



LOiLUk 



duj Lo^k-{k — 2(15) 
kik^LL!'^{k — LoY' 

kiijj'^{k — 



1 



k-iU^ 



12 



k- (k — Lu) 
a;2(fc-c3)2 



(fc — (15)2 — aj4(ti)4 — k/i) 



k^uj^{k~ujy 



fc2 

1- — 



k^uj^ik-csy 



LOiLOm 



duj tji{uj)tnk{k - LU) 



k^LU^{k-Uj)^ 

[{2uj - k)i5kj - 2kk6ij + 2kjSik] [{2uj - k)m5in + 2fci(5„ 



(2.26) 
(2.27) 



(2.28) 



k-{k — w) 

L3^{k-LUy 



(2.29) 



(2.30) 



In the equation, (j2.22p . for Taa (the gluon polarization), the first two integral terms (which contain the possible 
energy divergences at one-loop) cancel explicitly and this is due to the cancellation of the Faddeev-Popov determinant 
against the determinant arising from the Gaufi' law constraint that is the raison d'etre for the first order formalism 
used in this work fl']. 

Let us now discuss in more detail the equations (|2.24f[2730|) . Under the transformation uji — s- —UJ4, the right-hand 
side of Eq. (|2.24p and the second term of the right-hand side of Eq. (|2.28p change sign which means that these integrals 
must vanish: i.e.. 



^jili^l^ — r^|(fc4, P) — 



(2.31) 
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and 



P k-(k-Lu) 



(2.32) 



ki^k^uj^lk — 

Comparing equations (|2.25p and (|2.26p . we have clearly that 

rW(fciP) = rW(fc2) (2.33) 

which is also a consequence of the cancellation of the Faddeev-Popov determinant. The tensor equations, (|2.27p and 

p.30|) . must be decomposed to extract the dressing functions. We notice however that contracting Eq. (|2.27p with 
Ikiik) gives us the more useful combination 

ri'M,k') + Ti%lP) -r(i)(P). (2.34) 

On the other hand, contracting Eq. (|2.27|) with tki{k) gives 



(d-i)rW(fc|,fc2) = iv,^^ 

Turning to Eq. (j2.30p . contracting with lmi{k) gives us 





k-w^' 


1 u;2(fe_j)2 


k'^uj^ 



k^uj^ik^ujy 

whereas contraction with t„ii{k) gives 

(d-i)rW(fc2^fc2) ^ 



doj k-{k — 2uj) 

fc4[j2(fc-Cj)2 



^ ^ fc2 k.{k-Lo) 



L0^{k-UJY 



(2.35) 



(2.36) 



/A, Nc ^ 

k^U}^{k — LLlY ^ 



duj 



k^LLi^{k — a;)2 

[(2a; - k),6kj - 2kk6,j + 2kjdik] [{2uj - k)mSin + '^h^mn - '^knSmi 



t„ii{k)tji{uj)tnkik -LO) X 



(2.37) 



Expanding the transverse projectors of Eq. (|2.37p . contracting indices and where possible canceling denominator 
factors, we get 



(d-l)r«(fc42,fc2) ^ Ncfl' 

1 



duj 



k'^Lu'^{k — (15)2 
duj 

Puj'^ik-Luf 



1^4 lr2 



{d-2)L0^ + 



/fc2 



(6 - 4d)fc2 + 2{d - 1)-^ + 2(1 - d)oj^ 
k^ 



-k'^ - -k^k-u + Ud - 5)k-uj'^ - 2^- 
4 2 fc2 



1 



fc6 



8a;2(fc- 



(2.38) 



In order to reduce the number of integrals that we need to compute, we use the identities 



k-uj = l/2{k^ +uj^ - {k~ujY)-k4L04, 

, ,2 



2 ^2 
UJ — UJ 



to get 



(d-i)rW(fc2,fc2) ^ 



Nr. 



duj k* 



uj'^ik — u!)^uj'^{k — a;)2 

(4rf_5):^_fc2-fc2_^ 

/c2 fc4 



fc2 fc"' 

l + (10-8d)— + 3— 

fc2 fc4. 



duj 



Nc , 



iu;2(A; — a;)2a;2 



a;2(fc — a;)2tJ2 



+Ncii' 
+Ncti' 



duj 



uj'^{k — Lo)'^ 
duj 

uj'^{k — u;)2 



k 



.k- 



1 + (3 - M)— + 2— + 2— _ - 1 + 2(1 - d) — 



/c2 k^ 



fc2 Vfc2 



/c2 fc4 



3d- 



17\ fc2 5fc'' fc-cJ /7 



4 / fc2 4 fc4 /,2 



3fc2 



2 2A:2 



+ (rf-2) 



to 



A:2 
(2.39) 
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We must also deal with the integrals for the T^^, Tao- and Taa Green's functions, given by equations (|2.29p . (|2.32p 
and (|2.36p . It turns out that all these expressions are related as we will now show. Starting with TAa, Eq- (|2.32p . we 
can easily show that 



r^Aliklk')^id-^)N,^l' 



duj uJik-{k — 2uj) 



duj (2a>4 — k4)k-Lu'^ 
k,^k'^uj'^{k — uj)'^uj'^ 



(2.40) 



The second integral expression will be reduced below but it is convenient for now to leave it in its present form. Now 
consider Eq. (j2.29p . If we first recall that the last factor was originally written in symmetric form 



then it is easy to show that 



tij{uo)tji{k - uj) 



duj (fc4 — 2^4) 

kik'^uj'^{k — LoY 



^ ^ k^ ^ k-{k-Lu) 



{d~l)LO^-e 



duj (2CJ4 — fc4)fc-cJ2 



ki^k'^Lo'^ik — Lo)'^ijj'^ 



(2.41) 



(2.42) 



We immediately recognize the latter factor as occurring in Eq. (j2.40p . As for the former, firstly we have as a general 
result that 



duj (ki — 2CJ4) 



= 0. 



LL!^{k — Uj)'^ 

It can also be shown (using the results of Appendix [A|) that 

duj (fc4 — 2^4)0}'^ f duj uj/^k- [k — 2lo) 



k^k'^bj'^ik — LoY 



k4k^uj'^{k — toy 



Thus we see that 



rW(fc|,P) = r«(fcI,^^^)- 



Turning now to Eq. (|2.36p . it is easy to show that 



T'-AAiklP)=-N,fi' 



and further manipulating, we get that 



(d-l) f duo k-{k-2ujf 



k^uj^ik-ojf 



duj ^4(20^4 — k4)k-u}'^ 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



Taa (^4 7 ^^) 



k^ 



^ ^Aai^l, k^) 



(2.47) 



Let us finally reduce the expression for Taui Eq. (|2.40p . to the set of most basic integrals as we did for V aa- Using 
the same techniques as before, we have that 



^'^Aaik'ii k^) 



1 



duj 



dtU k4(jJ4 



1 k* 



U!^(k — LoY'LJ 



duj 

U^k-Lu)'^ 



d - 1)-^ - 1 + 2— - 2 d - 1 ^ + 2— 



1p 2 p 



(2.48) 



Before proceeding, let us briefly summarize the results of this section. We have written the various one-loop, two- 
point proper Green's functions in terms of the most simple collection of integrals. and T^aL trivial, Eq. (j2.3ip . 



We must calculate T^P, Eq. (|2:25| . which will also yield V^^l via Eq. ((233| . With T^^ calculated using Eq. ((235)) . 



(1) 



-(1) 



we also get T^^l from Eq. (|2.34p . FJ^^ is calculated using Eq. (|2.48p which then gives us F^V and T'^aa from relations 
(|2.45p and (|2.47p . respectively. Finally, we have to calculate F^]^ using Eq. (|2.39p . The necessary integrals will be 
derived in the next section. 



-(1) 
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3. NON-COVARIANT LOOP INTEGRALS 

In this section we consider the nontrivial loop integrals that arise in this study (i.e., massless one-loop two-point 
integrals). Quite generally, the integrals that arise within Coulomb gauge perturbation theory can be classified into 
two categories - those that can be evaluated using standard techniques such as Feynman parametrization or Schwinger 
parameters (and which are detailed in Appendix |X| and those that cannot. The latter category clearly requires a 
different approach and to this effect we derive a technique based on differential equations and integration by parts 
[IBP] suitable for the non-covariant setting here. We consider the three integrals: 

A{kl,k'^) ~ / ~277 \2 ^2 ' ('^■^) 

,2 r2\ f 



B{ki,k') = / ^ . (3.3) 

J LU^{k~CjyUJ^{k~Uj)^ 

It is convenient to introduce the following notation: x = k\, y = k"^ , z = x/y, v — y/x. We find that the above 
integrals can be written in the form 



jx + y)-^-'' 
(An) 



^i^^y) = Z'Le /"(^)' (3-4) 



A^ix,y)^kUix,y) = k^-^±^I^Uiz), (3.5) 

Bi-,y) - ^"(t)t-r' A(^) (3.6) 

where the fi{z) are functions of z and which may contain an ultraviolet divergence in the form of a simple pole as 
£ ^ 0. The functions fi{z) can (and will) be written down in analytic form for e — > 0, but it turns out to be more 
useful to write some of the various parts in terms of an integral representation and also as asymptotic series in z or 
V which are more amenable to eventual numerical evaluation. 

A. Derivation of the Differential Equations 

To begin, let us derive the differential equations that these integrals obey. Consider the general integral (n = 0, 1) 

Since /" is a function of two variables, there are two first derivatives: 

= /^^^f.^iai^j, ,3.8, 

OK4 J Lo^ik — ujYuj-^ y [k — LoY J 

5/" _ r duo Lol { k-{k-uj)' 
'Wk ~ J uj^{k-uYQ^ 1 [k-uf 



Now, there are also two integration by parts identities: 

0=/-^.^^ ,. - / ..""l.J - + l-24-2^^t^L (3.10) 



dwi Lu'^{k ~ LoYuj^ J LL!^{k — ujYuj^ \ LU^ {k — Uj) 

n- f rt ^ ^'^4 _ f duuj'i ^ ^uj^ ^ uj-{uj-k) \ 

^~ J du.iu^k-ioYco^ ~ J u^k-coYcu^Y {k-LoY j' 

Adding these two expressions gives 

^ f duj ( , „fc-(fc — w)l „x 
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from which we have the important identity 



Qjn gjn 

fc4Tn- + = (rf + n - 5)/". 



Expanding the numerator factor, Eq. (|3.12p can be rewritten 



= 



d + n-A 



Similarly, Eq. p.lOp becomes 







We can now rewrite Eq. (|3.8p as 



dr 

^ dkA 



n- 1+2— 



^4 



and so we arrive at the temporal differential equations for A and A"^: 



hi 

^4 



dA 

dki 



dkA 



l + 2{d-4) 



2{d-3)§ 



^4 



A + 2 



duj 



u;^{k-ojf 

duj UJ4 



- 2 



duj 



- 2 



LJ^k^Uj)^' 

duj uji 



(3.13) 
(3.14) 

(3.15) 

(3.16) 

(3.17) 
(3.18) 



The differential equations involving the spatial components are subsequently given using Eq. (|3.13p . 
Now let us consider integrals of the form B. Again, we have two IBP identities: 



0= / d 



d 



du 



diUi uj^(k-uj)^uj^{k-2ujy 

_d uji 

duJt uj'^{k-uj fuj'^{k-2ojY 



duj 



uo^{k-uYuj^{k-2ujY 

duo 



w| _ ^ 0^4(^4 - fc4) 

u;2 {k-uf 



(3.19) 



d-2-2— -2 



(w — fc) — k') 



(k-ujy (fc-a;)2 



and by adding the two we see that 

: 



duj 



uj^{k-ujYuj^{k-2ujY 
from which we have the important identity 



„k-(k — uj) ^k-(k — uj) 
d-7 + 2^ r:f + 2^ 



dB OB 
'^4Tn- + fcfc7Ti- = {d-7)B. 
oki okk 



Rewriting the definition of B as 

2S = 2 

and using Eq. (|3.19p gives 
B-2 



Cj4(fc-tj)2(I;2(fc-2a;)2 



■I 



duj 



du 



duj 



oj^{k~ojY{k~uY 

LUi{uJi — k^) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



w2(A;-w)4j2 7 ,^2(fc_^)2^2(fc_2tJ)2 [ {k - Uj)^ 



-2- 
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However, with Eq. (|3.14p we obtain the temporal differential equation for B: 



dkA k"^ k"^ 



duj 



(3.25) 



and the spatial equation is given via Eq. (|3.22p . 

The differential equations are best written by evaluating the standard integrals in terms of e (see Appendix [A| and 
with the notation x and y for the temporal and spatial momentum components. The result is [A^ = k^A): 



2x 
2x 



OA 

dx 
dA 
dx 



1 - (2 + 4e) 



l - 4e- 



x + y 



A~2- 



X + y 



X + 4:{x + y) 



-l-e 



x + y 



A-2- 



x + y 



X + 2ix + yy 



2x— = -'iB + U + ^e){x + y)-^A- A- 



where 



X = 



The spatial differential equations subsequently read: 



Y = 



(47r 



'r(i-2£)- 



2?/ 



dA 
dy 



„ dA 
dy 

„ dB 
dy 



-3 - 2e + (2 + 4e) 



-1 - 2e + 4£- 



x + y 



A + 2 



y~ 



x + y 



A + 2-^ 



x + y 



-X~4ix + yy 

x + y 

X-2{x + y)-'-'Y, 



-l-E 



Y. 



= (-1-2£)B- (4 + 8e)(x + j/)-iA + 4 



-l-e 



X + y 



-X. 



(3.26) 
(3.27) 
(3.28) 

(3.29) 

(3.30) 
(3.31) 
(3.32) 



B. Solving the Differential Equations 



The differential equations 
solution. Taking A to start (the other two follow the same pattern), we can write 

A{x, y) = Fa{x, y)GA{x, y) 



have a rather special structure that allows for a relatively simple 

(3.33) 



such that 



2x 
22/ 

FA{x,y)2x 
FA{x,y)2y 



dFA(x,y) 

dx 
dFA{x,y) 

dy 
dGAix,y) 

dx 
dGA{x,y) 
dy 



1 - (2 + Ae) 



x + y 



^3 - 2e + (2 + Ae) 



FA{x,y), 

FA{x,y), 



x + y 



x + y 



X + A{x + y) 



-l-s 



= 2 



x + y 



X - 4:{x + yy^^'Y. 



The solution for the first pair of equations can be written down without difficulty and is 



(3.34) 
(3.35) 
(3.36) 
(3.37) 

(3.38) 



By inspection, the constant C = since the function FA{x,y) has the dimension [y] ^ (i.e., the function Fa carries 
the dimension of the integral A). The remaining differential equations for Ga now read 



1/2, -l/2+s(^ + y)-l-2^2x^^^^ 



X ' y 



dx 



x + y 



X + 4{x + yy 



'Y, 



l/2„-l/2+e(^ + y)-l-2^22/^^^i^^ 



dy 



x + y 



X~Aix + y) 



Y. 



(3.39) 
(3.40) 
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Because the derivatives with respect to x and y are distinguished only by a minus sign, the function Ga rnust be a 
dimensionless function of the dimensionless ratio v — y/x or (equivalently) z = x/y. We choose to express Ga as a 
function of v to avoid singularities and the two partial differential equations collapse into a single first order ordinary 
differential equation: 



dGAjv) 
dv 



Using the integral representation of the hypergeometric function [T^: 



V 



gives the solution 



Ga{v)~Ga{Q) 



1 



1 

h\l + vY' 

^l/2-2e 



2F1 a,l;l + 6; 



(3.41) 



(3.42) 



l/2-2£(l 



\-2e- 



X2F1 -2e,l;3/2-2e; 



— — -s, 1; 3/2 - e; -— 
1 + \ l + v J 



(3.43) 



(3.44) 



1/2 -e(] 

Let us now show that the constant Ga{0) = 0. The general solution for A(x, y) as ?/ — ^ would have a term 

Fa{x, y ^ 0)Ga{0) ^ lim x^^^y-'/^+' 

but in the original partial differential equation p. 261) . which must still be defined for all values of y (although the 
coefficients may be singular) there is no such term and the constant G'yi(O) must vanish. We thus have the solution 

A{x,y) ^ 



{x + yy 



X 



1/2- 2e Vl 



2F1 -2£,l;3/2-2e; 



1 



2Y 



1/2- e' 



2F1 -£,l;3/2-e; 




(3.45) 



As it stands, the solution is not of much use since the hypergeometric functions arc somewhat cumbersome. However, 
we are primarily interested in the solution as e — > 0. In this case 



2F1 ~e, 1-3/2- e; 



= l + ef{v) + Oie^) 



and upon expanding A in powers of e the functions f(v) happen to cancel such that 



A{x,y) 



(x + y) 



'l-e 



(47r 



\2-e 



7 



41n2 + 21n 



x + y 

y 



0{e) 



(3.46) 



(3.47) 



A few remarks are in order here. The overall dimension of the integral A can be written in terms of the covariant 
factor x-\-y as for the integrals in linear covariant gauges and result in the standard logarithmic factor, singular on the 
light-cone and with a branch-cut extending into the timelike region. The ultraviolet divergence characterized by the 
1/e term is similarly covariant. The noncovariant component is logarithmically singular at y = 0. The cc-dependence 
of the noncovariant component, however, has the logarithmic singularity at a; -I- y = 0. This is actually rather crucial 
since otherwise it would be difficult to justify the analytic continuation between Euclidean and Minkowski space 
without further cancellations in forming the two-point functions. Lastly, it is remarkable that the integral A which 
defies standard evaluation techniques reduces merely to a combination of logarithms. Unfortunately, this simplicity 
will not be present in the integrals and B. 

Let us now turn to the function A governed by equations (I3.27P and (|3.3ip . Writing A{x^y) — F-jG^ as before 
leads to four partial differential equations: 



2x 



2?/ 

F^{x,y)2x 



dF-j{x, y) 

dx 
dF-j{x,y) 

dy 



1 -4e- 



x + y 
-1 - 2£ + 4£ 



FAix,y), 



X 



dG^jx, y) 



x + y 



-X 



x + y 
2{x + y)-^~'Y., 



FAi^,y)^y 



dG^jx, y) 
dy 



x + y 



X-2{x + y)- 



(3.48) 
(3.49) 
(3.50) 
(3.51) 
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The solution to the first pair is 

F^{x, y) = x-^'^y-^'^+'{x + y)-^', (3.52) 
the constant vanishing on dimensional grounds as previously. The latter two differential equations now read 

X-V2y-l/2+e(a; + y)-2e2a;^ 



dG-7-(x,y) y , , i ^ 



dx 



x + y 



X + 2{x + y)-'-'Y, 



^-l/2y-l/2+e(^^y)-2e2y 



dG^{x,y) _ „ y- 



dy 



x + y 



X - 2{x + y)-^-'Y. 



(3.53) 

(3.54) 



Again, as before the relative minus sign means that is a function of the ratio, but this time it is useful to use both 
V and z, the reason becoming clear shortly. In terms of w, we have that 



dv 



whose solution is 



X 



^;l/2-2e 



= t;-V2-2e(l + ^)-l+2ej^ _ „-l/2-e(i + „)-l+ey 



^1/2-e 



(3.55) 



1/2 -2e (l+vy 
whereas in terms of z we get that 

dG^iz) 



^,F. (1 - 2e, 1; 3/2 - 2.; ^ J - ^ (TT^^^^ ' ^' ' 



dz 



1 + 

(3.56) 



(3.57) 



whose solution is 



,1/2 



1/2 -2£ (1 + ^)1-2^' 



2F1 l-2£,l;3/2;■ 



+ 



2y 



1 + zJ l/2-e(l + ;j) 



l^.F.(l-M;3/2;^) 



(3.58) 



The general constant terms in both solutions vanish as before. We can now write the result for the full integral 
A'^{x, y) = k4,A in two ways: 



A\x,y) 



kA 



(x + y) 



l+e ■ 



X 



1/2 -2e \ l + v 



Fi l-2£,l;3/2-2e;- 



Y 



5^1 (^l-£,l;3/2-£;^^|, (3.59) 



A^{x,y) 



{x + yy+^ 



1 + vJ 1/2-e 

|-2X(1 + zr2Fi (1 - 2e, 1; 3/2; + 2Y2F, (l - e, 1; 3/2; | . (3.60) 



Now, in contradistinction to the integral A, the expansion of the hypergeometric functions in powers of e does not 
yield a simple result - quite the contrary. However, one can expand in powers of either v or z and then take the limit 
e ^ which will prove useful in later numerical evaluation. Actually, in the case of the expansion with v, one must 
first collect together terms involving lni;/(l + v) since these factors do not have an expansion around v = 0. The 
results are ultraviolet finite and are: 



A\x,y) k,^^ + y^ 



-l-e 



A'^{x,y) ''=° ki 



(47r)2-£ 
'J - 

{x + y)-^-' 



In 



l + v 



21n2 



^ 4 404.44 

-2 v-\ V H V 

3 15 35 63 



20 124 , 

+4H V v'^ + 

9 225 3675 



3 3244 4 



(47r) 



2-£ 



In 2 



^+3^-15^ 



19845 



■}■ 



(3.61) 



35 



63 



3" -15" 



105 



23 
378' 



(3.62) 



The statement that the expansion of the hypergeometric functions in powers of e is nontrivial is fortunately not the 
whole story. It is possible to find the full solution as e ^ by rewriting the hypergeometric functions in their integral 
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representation, expanding the integrand and exploring whether or not it is possible to evaluate the resulting integrals. 
In the case of using v as the variable this is not the case - the expression reads: 



A^{X:V) = fc4 



[x + ij)-^-^ (1 + v) 



2-£ 



r dt ft 

-41n2arctan(-v/i;) — / In 

^ Jo Vtil+t) \l + t 



(3.63) 



and the integral cannot be done directly in terms of known functions (or to phrase it more properly, the direct result 
is not known to the authors at present although indirectly the result can be inferred from the following discussion). 
However, with z as the variable, we get 



(x + y)"^-" (1 + z) 



and the integral is 



dt 



Vtil + t) 



(4^)2-^ ^ 
ln(l + i) = 7rln2 - zln(yz- z) 



4 In 2 arctan ( \/z) — 



dt 



Vtil + t) 



In (1 + t) 



(3.64) 



In 2 + In (1 + z) - In (1 - iVz) - ^^t^ iV^ - ^) 



+iln iVz + %) 
'1 



In 2 + In (1 + z) - In (1 + %Vz) - -j^^ iV^ + «) 



idilog( - Vz] + idilog 



(3.65) 



in terms of the dilogarithmic function |Ta | . An expansion in powers of z yields the expansion Eq. (|3.62p which serves 
as a useful check. Whilst it is gratifying that the integral can eventually be written in terms of known functions, 
these functions must still be evaluated at some stage. For this reason, the integral form Eq. (|3.64[) and the asymptotic 
forms Eqs. (|3.61l3.62p are actually of more practical use. 

It is possible to discuss the analytic structure of A^. Clearly, it is ultraviolet finite. With the expansion Eq. (|3.6ip 
in V = y/x we see that as y ^ there is a logarithmic singularity as seen previously for the integral A. As for the 
singularities in x, by rewriting the integral 



dt 



V^ Jo Vtil + t] 



■ln(l + <) = 



dt 



Vtil + zt) 



In (1 + zt) 



(3.66) 



and knowing the analytic properties of arctan (^/z) / V^ '^^ that the singularity occurs for z — —1 (i.e., for x+y — 0) 
with branch cuts extending into the timelike region. Again, this behavior is the same as for A. 

Finally let us discuss the integral Bix, y) satisfying equations (|3.28p and p.32p . We separate the function into two 
parts as previously: 



Bix,y) = Fb ix, y)GB ix, y) 
such that Fb and Gb obey the following differential equations: 

dFBix,y) 



(3.67) 



2x 

2y 

FBix,y)2x 
FBix,y)2y 



dx 
dFBix,y) 

dy 

dGBix,y) 

dx 
dGBix,y) 
dy 



-3FBix,y), 


(3.68) 


-il + 2e)FBix,y), 


(3.69) 


(4 + 8e)(x + y)-^Aix, y) - 4^^X, 

x + y 


(3.70) 


-(4 + 8e)ix + y)-^Aix, y) + 4^-^X 

x + y 


(3.71) 



The solution to the first pair of differential equations is 

FBix^y)^x-^'^y-''^-^ 



(3.72) 



with the possible constant vanishing as before on dimensional grounds. The second pair of equations now reads 



dGBix,y) 



j;'V2y-i/2-e^^ ^^myyj ^ {4 + 8e)ix + y)-^Aix,y)^A 



dx 

-Z/2y-l/2~eny dGBix,y) 

dy 



^iA + %e)ix + y) ^Aix,y)+A 



y-l-e 

X 

x + y 

y-l-e 



x + y 



X. 



(3.73) 
(3.74) 
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We are now faced with a potential problem - for general e, the function A(x, y) is itself a combination of hypergeometric 
functions and we have little chance of solving the resulting differential equations. However, we can write down A{x, y) 
for vanishing e. Again, the derivatives with respect to x and y are distinguished by a minus sign and we can rewrite 
the two equations as a first order differential equation in z = x/y. In this case, we cannot use the variable v since 
this leads to non-integrable singularities. After expanding in powers of e, the equation is 



dGsiz) _ 2(1 + 2e) 



dz 



(47r)2-£ 
2 

"(47r)2-^' 



1/2 



(1 + ^)' 



-2 



.1/2 



- - 7 - 2 In 2 

e 



-41n2 + 41n(l + z) + 0(£) 
0{e\ 



(3.75) 



Since there is no interference between the integration over z and the expansion in e, we do not need to consider the 
terms of 0(e) further. We thus obtain 



Gb{z) 



(4^)2 



4z3/2 
(1 + ^) 



i - 7-21n2 ) + 8(l + 21n2) 



L(l + ^) 



arctan (\/z) 



dtVi 

(TTTp 



ln(l + <) + 0(e) 

(3.76) 



where the constant Gb(0) vanishes as before. Using integration by parts, we have that 

dt 



-^ln(l + <) = - 



(l + t)2 "-^ ■ - 2 
and so, the full solution can be written 
(a: + y)-=y-2 



Vi{l + t) 



In (1 + 1) 



(1 + ^) 



ln(l 



(1 + ^) 



arctan (-^/z) 



(3.77) 



B{x,y) = 



(47r 



|2-£ 



-4(1 



- - 7-21n2 

s 



z{l + z) 



In (1 + z) + 



dt 



z Jo Vt{l + zt) 



+ In (1 + 
ln(l + zt) + C'(e 



16 In 2 



z(l 



23/2 



rctan (Vz) 



(3.78) 



The factor y~^ is extracted globally for convenience as it happens that the integral B(x, y) is multiplied by when 
forming the gluon polarization. Without the factor y"^ , the expression is finite as y — > 0. The singularities in x occur 
for z = —1 as before, with branch cuts extending into the timelike region. An expansion around z = is possible and 
reads 



-ey-2 



.^0 [x + y) 



+ 



16 

y 



-4 
46 



28 



z^ + 



\l-z 



z^-z^ + z^ + 



Sin 2 



7 o 9 4- 
9 11 



202 . 91 4 

z-^ z^ 

27 11 



(3.79) 



Comparing Eq. (|3.78p with Eq. p.64p and using the known expansion for A/^i^x, y) for small i;, Eq. (I3.6ip . we can also 
expand B[x,y) and the result reads: 



B{x,y) "= 



[x + y) "y 



e„-2 



(47r) 
-4 I In 



2-£ 



-4 



1 



1 

e 

2 In 2 



64 
— ? 
9 



368 
'75^ 



(3.80) 



Having derived the integrals, it is pertinent to see if they can be checked using standard techniques. Since this is 
a somewhat technical exercise that does not add to the discussion here, we present the details in Appendix |B] It 
is also useful to plot the integrals and their asymptotic expansions. The functions fi{z), defined in Eqs. (|3.4ll3.6p . 
and with the factors proportional to (1/e — 7) removed are presented in Figures [TJO All functions are monotonically 
increasing with z. One can see that the asymptotic expansions do indeed represent the functions within their domains 
of applicability. The functions faiz) and fiiz) are logarithmically singular as y — > (2 — > cxd with x fixed), whereas 
fb{z) is finite. All functions exhibit singularities as z ^ — 1 (especially strong in the case of fb although this can 
be attributed to the different choice of prefactor in the full integral B). However most importantly, there are no 
singularities at z = (or equivalently x = for finite y), such that the continuation to spacelike Minkowski space is 
entirely justified. 
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0.001 0.01 0.1 1 10 100 1000 -1 -0.5 0.5 



z z 

FIG. 1: UV- finite part of the function fa{z). Left panel: Euclidean values of z. Right panel: continuation of fa{z) into the 
spacelike Minkowski region. 




FIG. 2: The function fi{z) and its asymptotic expansions for z ^ and v = 1/z ^ 0. Left panel: Euclidean values of z. Right 
panel: continuation of f4,{z) into the spacelike Minkowski region. 



4. ONE-LOOP TWO-POINT FUNCTIONS 



Having now evaluated all the integrals that occur (see the previous section and Appendix E|, we can now return 
to the expressions for the one-loop, proper two-point functions and write out our results. To do this, it is convenient 
to define two combinations of functions: 



f(z) = 4 In arctan ^/z — / —= ln(l 

Jo Vt{l + zt) 



zt), 



g{z) = 21n2 - ln(l + z). 
[These functions are actually variations of the finite parts of /4 and fa from before.] The results are: 



r«(j/) = 

(47r)2-e 

(47r)2-e 



(47r)2-6 
4 
3 



- - 7 - In - 



28 8, „ 

h - n 2 

9 3 



-0{e) 



1 

T -7- 



ln(^ 



52 
~9 



^ln2- 



0{e) 



0{e) 



^^Aai^^y) = ^''a^i^^y) ^ --'^AA{x,y) 



2-e 



1 

- - 7 - 



In 



x + y 



(4.1) 

(4.2) 
(4.3) 

(4.4) 

(4.5) 



i - 6(1 + z) + 3(1 + z)g{z) + 1(1 + z)(l + 3z)f{z) + 0{e) 
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fi,(z), numerical 

; ^<^^)^ z small 

1 —- f^(z), z large ; 
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FIG. 3: UV-finite part of the function /6(2:) and its asymptotic expansions for z ^ and u = 1/z 
values of 2. Right panel: continuation of jb{_z) into the spacelike Minkowski region. 



0. Left panel: Euclidean 



(47r 



2z 



(1 + z) + (1 + .) -- + -- + 



y + — (l + ^)-3z(l + z) + -5(^) 



1 1 11 



4z 



(4.6) 



(4.7) 



It is immediately apparent that for finite y, there arc no singularities in any of the above at x = (or equivalently, 
2 = 0) since /(O) — 2g(0) = 41n2 cancels the 1/z pole. This is to be expected since none of the individual integrals 
or their prefactors are singular at this point. 

Let us now construct the one-loop propagator dressing functions. Writing Daf) — o'^a^^ + g^D^'^^ and using 
Eqs. (|2.31l2.33l2.45p . Eq. (PTT|) becomes 
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^AAi^^y) - 
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The last equation of (|2.1ip becomes an identity by virtue of the relation Eq. (|2.47p . as it should. Putting in the above 
results for the one- loop, two-point proper functions, we have (in Euclidean space and in the limit e — > 0): 
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A few remarks are in order here. The momentum dependence of the relationship between the vector propagators and 
proper two-point functions, Eq. (|2.1ip . is such that the only occurrence of a momentum dependent UV-divergence is 

within the D-^^^ propagator ~ the factor 1/(1 + z) is otherwise canceled within the combination fJ^],-!- F^^i. Indeed, 
the divergence structure of the one- loop, proper two-point functions has been known for some time [l|. The momen- 
tum dependent coefficient of the 1 / e-pole in d'^^ is symptomatic of the fact that the Tr-field is not multiplicatively 
renormalizable since the 7r-field has its origins in the linearization of the (composite) chromoclectric field term of the 
action that is central to the first order formalism. However, we are further able to see that for the UV-finite parts, the 
kinematical singularities on the light-cone (z — —1) reside purely in the logarithmic term and the functions /(z) and 
g{z) which are logarithmic in character. There are no singularities in the Euclidean or spacelike Minkowski regions 
(z > —1). Hence, we can conclude that the analytic continuations between Euclidean and Minkowski space have 
entirely the same character as in linear covariant gauges - that is to say that the continuation is justified. 

Although it is not our intention to discuss the renormalization aspects of the two-point Green's functions in Coulomb 
gauge, at the one-loop perturbative level it is possible to identify two renormalization group invariant combinations 
of propagator dressing functions via the coefficients of the 1/e poles. In Landau gauge, a renormalization group 
invariant running coupling may be defined through the combination of gluon and ghost propagator dressing functions: 
q^DaaD^, This stems from the Slavnov-Taylor identity which expresses the universality of the coupling and 

the Landau gauge property that the ghost-gluon vertex is UV-finite. Given that at one-loop, the coefficient of the 
1/e pole of is the first coefficient of the /3-function (60 = — llA''c/3(47r)2~^) and is gauge invariant, the combination 
DaaD^ in Landau gauge has a 1/e pole with the coefficient llA'c/3(47r)2~^. In Coulomb gauge, the results above 
clearly show the same result. However, the individual coefficients for Daa and Dc are different from Landau gauge. 
The second renormalization group invariant combination of propagator dressing functions is particular to Coulomb 
gauge and is g^D^^, |l]. The coefficient of the 1/e pole in Dg-a above clearly confirms this. 



5. SUMMARY AND OUTLOOK 



A one-loop perturbative analysis of Coulomb gauge Yang-Mills theory within the first order formalism has been 
undertaken. The various propagator and two-point proper dressing functions have been explicitly evaluated at this 
order. In order to do this, dimensionally regularized results for the noncovariant two-point loop integrals inherent to 
Coulomb gauge have been derived using techniques based on differential equations and integration by parts identities. 

The results for the two-point functions are rather interesting. The dressing functions are dimensionless functions 
of two independent variables and a mass scale introduced via the regularization. These functions can be split into 
two parts - one unambiguously connected to the UV-divergence involving also the logarithmic behavior normally 
associated with covariant gauges and a second part which is a UV-finite function of the ratio of the temporal to 
spatial components of the momentum. 

The analytic continuations between Minkowski and Euclidean space within the noncovariant setting can be justified 
on the grounds that the possible singularities occur on the light-cone, with branch cuts extending into the timelike 
Minkowski region. It is seen explicitly that for spacelike Minkowski and Euclidean momenta there are no singularities 
in any of the two-point Green's functions, which is as it should be. 

The outlook for future work done in Coulomb gauge is rather promising. The most direct continuation of this work 
is to consider the vertex functions of the theory. A generalization of the differential equation technique to the various 
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one-loop three-point loop integrals certainly appears feasible albeit challenging. Subsequently, a two-loop perturbative 
analysis studying, for example, the cancellation of potentially energy divergent integrals or the renormalization would 
be of great interest. A second area of interest would be to include quarks and to study physical high energy processes, 
as has been done in linear covariant gauges. The relationship between the covariant and noncovariant descriptions of 
the same phenomena will undoubtedly lead to greater insight into the physical mechanisms at work. 

One of the motivations for studying Coulomb gauge is that nonperturbative phenomena such as confinement and 
bound states may be better understood in this gauge. The analysis of nonperturbative physics is however greatly 
constrained by the perturbative behavior. As an example, consider the evaluation of nonperturbative loop integrals - 
the phase space of the integration measure still contains the perturbative domain and renormalization is still necessary 
despite the fact that one may be considering infrared external momentum scales. Further, the techniques developed 
here to evaluate the perturbative integrals will almost certainly be of help when studying the Dyson-Schwinger 
equations nonperturbatively. Also, the perturbative expansion (although asymptotic) is of great use in verifying 
nonperturbative identities such as the Slavnov-Taylor identities. The Slavnov-Taylor identities are the focus of 
present work [l8l |. 
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APPENDIX A: STANDARD INTEGRALS 



There are two particular types of integral that we wish to consider in this appendix and both can be done using 
standard techniques. We use the Schwinger parametrization method [ll| . The integrals considered have two or three 
denominator factors (with arbitrary powers), at least one of which contains both spatial and temporal components. 
We list results for all the possible vector and tensor integrals that arise. 

Consider then the integral (in Euclidean space) 



Using the identity [16| 

1 1 



a'' r(i/) 

we have 



/"OO 

/ daa^-^expl-aa}, ^ly > 0, SRa > (A.2) 
Jo 



r(M)r(i^) Jo 

Shifting variables 



J = ^ / dad(3a''-^f3''-^ f dujexp{~{a + l3)w^ + 2(3k-uj-(3k'^}. (A.3) 

i (^)l (l^) Jq J 



uj^uj-\ (A. 4) 

a + p 



gives 



ri^i)r{iy)Jo ' J V ^^13 

and rescaling, lj ^ {a + P)~^^^uj, leads us to 



/ = ^ . .^^ . . / dad(3 af"-^ p''-^ f duj exp\~{a + p)w^ - -^kH (A.5) 



r(/i)r(z.) Jo 

The integral over lo can now be carried out 



^— - / dadpa''-'^(3''~\a + l3y~^expl ^k^\ f duj exp \-wH . (A.6) 

)r(z^) Jo [ a + l3 } J 
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By inserting the identity 1 = dXS{X — a — f3) and rescaling a Aa, /3 ^ A/3 we have 

^ ^ l'\adpS{l-a-p)a'''^p''-\a + (3Y~^ /°° dA A^+"+^-^ exp ( -A^^fc^ 



(4^)2- r(^)rH ^ ^ ' ^' ' " ' Jo n « + 

— - / daa''-\l-ay-^ dXV'+''+'~^ exp {-Xa{l - a)k^} . 

')r(i^) Jo Jo 



(47r)2-£ r(^)r(!. 

The integral over A can be done and is a variation of Eq. (jA.2p giving 



(A.8) 



(47r)2-- r(M)r(i/) Jo 

Finally the integral over a can be done (it has the integral form of the beta- function) to give the final result: 

duj _ [P] '^^"^'^^ TQi + iy + e~2)T{2-fi- e)r{2 ~ v - e) 

Now let us consider the vector integral 

/duj uji 

In order to proceed we notice the following. Under the change of variables uj uj + k(3 / {a + P) we have: 
duj exp{-{a + (3)w^ + 2pk-Lu - (3k^} ^ J duj exp ^^-{a + (3)w^ ^'^^}' 



2(3 i dujiuj4-ki)exp{-{a + P)w^ +2pk-Lj^ (3k^} =-2l3 I dij —^k^exp I ^{a + P)w^ ^ -^^k'^ \ (A.13) 



Differentiating with respect to ki gives 

i^2nk-i,, ^ Ftk^\ = -9.Ft I d„, , 

a + P I a + /3 

which shows us that 

ffw W4exp{-(a + /3)w2 + 2/3/c-a;-/3fc2} = f dcj —^k^ex-p l-{a + P)w^ - -^^kH. (A.14) 

Further differentiation gives rise to expressions for integrals involving other numerator structures. Proceeding as 
before, we have the results 
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Using exactly the same techniques we have 
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APPENDIX B: CHECKING THE NONSTANDARD INTEGRALS 

Since the integrals A, and B must be derived using nonstandard techniques, it is worthwhile checking them 
where possible against available results. It turns out that the expansions around z = Q, Eqs. (j3.47|3.62|3.79p . may be 
checked analytically. An expansion around w = is not possible, since all integrals are divergent as y ^ 0. 

To begin, consider the integral A, Eq. (jS.ip . Using Schwinger parameters [Tlj , we can rewrite the denominator 
factors as exponentials, the result being: 
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Scaling the integration variables aj4 ^ (a + /3) ^/2w4, tJ ^ (a + /3 + 7) ^/■^o; then allows us to do the momentum 
integration, leaving the parametric integral 
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This last equation we denote as the parametric form of the integral. For general values of z, it cannot be done because 
of the highly nontrivial denominator factor (and clearly which is why the differential technique has been developed). 
However, knowing that the result is well-defined at z = 0, we are able to make an expansion around this point and 
then do the resulting parametric integrals. To second order in powers of z we have: 
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The parametric integrals can be done without difficulty and the result is 
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which agrees explicitly with Eq. p.47p . Actually, with hindsight and given patience in expanding the integral, it 
appears that this particular integral would be possible by resumming the series expansion! 

Turning to the integral A"^, Eq. (|3.2p . the parametric form of the integral can be written down almost immediately 
given the previous case. It reads: 



Expanding again to second order in z gives 
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which agrees explicitly with Eq. (|3.62p . Although at this order it may appear reasonable to suppose that one may 
resum the series to recover the full function, a quick glance at the full expansion Eq. p.62p and the solution Eq. (|3.64D 
tells us otherwise. 

As might be expected, the integral B, Eq. (|3.3[) . is rather more complicated. The parametric form reads: 
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The last factor can be expanded in powers of z and to second order, we get 
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The integral over P can be done and gives 
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Now, in order to do the last pair of parametric integrals, we change variables with a — r(l — s), 7 = rs such that now 
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The two-dimensional integral in a and 7 is now separated into two parts which can be done in turn. The integral 
over s yields powers of r and (1 — r) and the integral over r subsequently leads to the familiar combinations of gamma 
functions. Expanding in e and completing the expansion of z by including the prefactor 1/(1 + z)^ we finally arrive at 
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which agrees explicitly with Eq. (|3.79p . 
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